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£N) Abstract. We report on multifrequency radio observations 
CN made with the Effelsberg 100-m telescope, the IRAM 30-m 
telescope and the Green Bank Interferometer between 1992 
^SJ and 1994 of the 7-ray quasar 0528+134. We present a new 
OO VLBI based map of 05 1 8+ 1 34 at 22 GHz with sub-mas angular 
^ resolution observed in November 1992. At that time the source 
was in a phase of brightening at all of our observing frequen- 
ts-} cies above 3 GHz. The increase of brightness may be related to 
Q\ activity in the unresolved core component of the quasar. 
*^ The VLBI map at 22 GHz (epoch 1992.85) shows a one- 
Q pided core jet structure of ~ 5 mas length. A new component 
1 close to the core indicates an apparent transverse velocity of 
Papp <20 (Ho=100 km sec -1 Mpc -1 , qo=0.5). A second com- 
^ ponent, seen also previously at 8.4 GHz & 22 GHz, shows ap- 
C3 parent superluminal motion with /3 app =4A±U. Superluminal 
motion is expected since strong Doppler boosting in this source 
is required in view of the 7-ray luminosity and the variability 
timescale, which violate the compactness limit. 

In summer 1993 we observed a major outburst at all fre- 
quencies higher than a few GHz, which peaked a few months 
after a strong outburst in high energy 7-rays and showed the 
canonical time evolution of a cooling and/or expanding electron 
distribution. Our data indicate that the outburst in the EGRET 
range originated very close to the central object of the AGN and 
that a remnant of this outburst moved further outward in the jet 
until it became optically thin at radio frequencies after a few 
months. 

During the flare in July 1993 we observed with the Effels- 
berg 100-m telescope an unusually strong decrease of the flux 
density by about 50% at 4.75 GHz and 10.55 GHz and slightly 
less at 2.695 GHz. This behaviour is also seen in the monitoring 



data at 2.25 GHz and 8.3 GHz taken with the Green Bank Inter- 
ferometer (NRL-GBI). The event lasted less than three days at 
the higher frequencies and more than two weeks at 2.25 GHz. 
For the case that this event is related to the intense radio flare 
some geometrical effects like a small variation of the viewing 
angle of the quasars jet orientated very close to the line of sight 
are considered, but found to be an unlikely explanation for the 
observed behaviour. Alternatively, an extreme scattering event 
by a small dense plasma cloud in the line of sight is able to 
match the observed time lag in the lightcurves if we take into 
account the mas-structure of the source and different spectra 
of the components on the basis of their brightness in the VLBI 
maps. 

The importance of interstellar scattering is stressed as 
0528+134 is seen in the direction of the dark cloud Barnard 
30 located at 400 pc distance in the Orion complex and that we 
note frequent fluctuations of the order of 25% at 2.3 GHz / 2.695 
GHz, our lowest monitoring frequencies, while the variations 
at higher frequencies are more smooth. 

It is shown that the appearance of the extreme scattering 
event in the light curves of 0528+134 depends strongly on its 
VLBI structure and on the observed spectral appearance in our 
two frequency VLBI data. Due to these strong constraints our 
successful modelling provides the best available evidence that 
unusual variability behaviour of AGN may be caused by ex- 
treme scattering events and that localised (r< 1 AU) and dense 
(n e =100-1000 cm -3 ) plasma structures do exist in the interstel- 
lar medium of our Galaxy. 

Key words: radio sources:0528+134 - variability - VLBI - 
ISM: structure 
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1. Introduction 

The quasar 0528+134 is one of the strongest AGN in the hard 
7-ray range as observed by EGRET on board of CGRO (Hunter 
et al. 1993). It is a m = 20.0 red object located in the direction 
of the dark cloud Barnard 30 at the periphery of the A-Orionis 
complex which is at a distance of about 400 pc (Maddalena et 
al. 1986). The redshift of 0528+14 has been recently measured 
to be z = 2.07 by Hunter et al. (1993), which makes the source 
one of the most luminous objects in the universe. 

0528+134 is a known variable radio source and has been 
regularly monitored for about 20 years (e.g. Aller et al., 1985). 
After the launch of CGRO in 1991 it was one of the first new 
detection's by the hard 7-ray instrument EGRET. The source 
is located in the galactic anticenter close to the well known 
7-ray sources Geminga and the Crab. Although this area was 
also deeply observed by COS B, 0528+134 was never noted 
during that mission lasting from 1975 until 1982. Since the 
fall of 1991 regular observations with the Effelsberg 100-m 
telescope have been made in the frequency range 2.695 GHz 
to 32 GHz. The source has also been monitored by the IRAM 
30-m telescope around 90 GHz and less frequently at higher 
frequencies. Shortly after its detection in the 7-ray range the 
source increased substantially in intensity first at 90 GHz and 
32 GHz and with some delay also at lower frequencies. Reich 
et al. (1993) have already presented these observations until the 
end of 1992. However, the limited time sampling so far did not 
allow us to establish a clear correlation between 7-ray and radio 
bands. We discuss the radio lightcurves up to the beginning of 
1994 in section 2. 

Many of the 7-ray bright AGN are identified with super- 
luminal sources. 0528+134 was not known as a superluminal 
source due to an insufficient number of VLBI maps. VLBI ob- 
servations at 8.4 GHz and 22 GHz performed in 1990/1991 
(Zhang et al. 1994, hereafter referred to as Z94) showed a com- 
pact source which was interpreted by two alternative models 
either as a one or two sided core jet structure. In section 3 we 
present new VLBI observations at 8.4 GHz and 22 GHz done 
in 1991/1992, which now allow to discriminate between the 
different signatures of the VLBI components and to obtain a 
better understanding of the kinematics of the mas-structure. 

Section 4 comprises a general discussion of the results. We 
found a rather unusual strong effect in the radio light curves in 
July 1993 when the source flared at 32 GHz and 86 GHz. This 
event is discussed in some detail in section 4. There we inves- 
tigate the possibility of an intrinsic effect within the source to 
explain our observations. Alternatively we consider a scattering 
effect caused by the galactic interstellar medium in the line of 
sight to explain the observed radio light curves. 

2. Multifrequency radio observations 

Reich et al. (1993) have already described the multifrequency 
observing method used at the Effelsberg 100-m telescope to 
monitor variable sources as detected by EGRET. The observa- 
tions result in quasisimultaneous flux density measurements at 



2.695 GHz, 4.75 GHz, 10.55 GHz and - if weather allows - in 
addition at 32 GHz. The angular resolution of the telescope is 
4.3', 2.4', 1.2' and 0.45' respectively. Orthogonal cross-scans 
are made centered on the source and a Gaussian fit is applied 
to the averaged forward and backward scan taken in each di- 
rection. Small pointing offsets are corrected for by the analysis 
software when calculating the peak flux density of the source. 
A detailed error analysis from the Gaussian fit parameters is 
made and the scatter in the data obtained from the calibration 
sources during one observing run is added for the total error 
of the flux densities quoted. As listed in Table 1 the errors at 
2.695 GHz and 4.75 GHz are quite small typical around 0.1 
Jy, while they increase at higher frequencies. The observations 
for one source at all frequencies are typically completed within 
30 minutes. 3C286 serves as the main flux density calibrator at 
all frequencies assuming 10.4 Jy, 7.5 Jy, 4.5 Jy and 2.1 Jy at 
2.695 GHz, 4.75 GHz, 10.55 GHz and 32 GHz. In the case of 
0528+134 the nearby secondary calibrator 3C138, which is just 
4 degrees apart, is always observed before or after measuring 
0528+134. This gives an additional check on the accuracy of 
all high frequency observations of 0528+134, which in general 
may suffer occasionally by small opaque clouds in the line of 
sight. 

We have also taken data at 2.25 and 8.3 GHz from the 
NRL-GBI monitoring program to complement our data base. 
The purpose of this observing program is to study extreme 
scattering events as described in Fiedler et al. (1994). 

The method of observing with the Green Bank Interferome- 
ter was already described by Fiedler et al. ( 1 987b) and Waltman 
et al. (1991). The observing frequencies are 2.25 and 8.3 GHz, 
and the maximum baseline of the interferometer is 2.4 km. The 
source was observed on alternate days near the meridian. Cal- 
ibration was based on four sources but ultimately referred to 
3C286. The rms-errors of the derived flux densities are around 
2% at 2.25 GHz and roughly 5% at 8.3 GHz. 

Steppe et al. (1988) have already described the observing 
and calibration method at the IRAM 30-m telescope to monitor 
variable sources. For 0528+134 mainly 86 GHz observations 
have been made in 1993. Flux density errors are below 8%. 
Z94 presented closely sampled observations at 22 GHz and 37 
GHz made between mid 1991 until 1993 which complement 
the Effelsberg and IRAM data. The spectra at all wavelengths 
below 32 GHz are inverted. The lightcurves based on Effelsberg 
and IRAM data from the end of 1991 until mid the spring in 
1994 are shown in Fig. 1 . In Fig. 2a we give a detailed view of the 
Effelsberg data and compare them with the NRL-GBI data in 
Fig. 2b, which we show as averages over two weeks to improve 
the readability of the plot. Here the uncertainty in the Green 
bank data is not due to observational errors, but to day-to-day 
fluctuations of the source which may sum up to two or three 
percent. Unfortunately, for a period of two weeks in the first 
half of July 1993 no data at 2.25 GHz from the Green Bank 
Interferometer are available due to hardware malfunction. 

Fluctuations at 2.695 GHz of about 20% - much larger than 
the errors in the observed flux densities - are rather frequent 
during the time of monitoring and these fluctuations are sig- 
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nificantly larger than those at 4.75 GHz, 8.3 GHz, and 10.55 
GHz. There is an excellent agreement between the Effelsberg 
2.695 GHz and the Green Bank 2.25 GHz data. This is an ad- 
ditional indication that most of the fluctuations seen at this low 
frequency have to be considered as real variations seen on time 
scales of months. This is a strong indication that stochastic in- 
terstellar scattering is important in the direction of 0528+134, 
which is not unexpected considering the structure of the inter- 
stellar medium in the Orion region. 

For a few months in fall 1992, when we made our 22 GHz 
VLBI observations, the mm-flux densities are around 4.5 Jy, 
which is lower than the peak values seen in mid 1992 which 
range between 5.5 Jy and 7 Jy. At the time of the 7— ray flare 
in March 1993 the 32 GHz flux was at a high level of about 7 
Jy. Unfortunately 86 GHz data are missing during that period. 
After a short time of slightly lower flux density level at 32 GHz 
we note in mid 1993 the absolute flux density maximum rang- 
ing between 8.5 Jy and 9 Jy at 86 GHz and 32 GHz. At the 
time of highest flux density at 32 GHz the 86 GHz flux densi- 
ties declined already and this may be indicative of a time lag. 
The 10.55 GHz and 8.3 GHz flux densities remained roughly 
constant at a high level until the end of 1993. The evolution of 
this source during the outburst follows the canonical behaviour 
of appearing first at high frequencies, and then with decreasing 
optical depth being affected at lower and lower frequencies. 

A very unusual behaviour, however, is the dramatic drop 
in intensity at 10.55 GHz and at all of our lower monitoring 
frequencies at that time in July 1993. In Table 1 we have listed 
the detailed results of our observations made in 1993. As evident 
from Table 1 and Fig. 2a there has been a dramatic decrease 
in flux density by about 50% at 10.55 GHz and 4.75 GHz 
and slightly less at 2.695 GHz. This happens within four weeks 
although details are missing by insufficient sampling. Therefore 
the observed minimum in particular at 1 0.55 GHz might be even 
below what we have noted for July 17. The rapid increase of the 
10.55 GHz flux density between July 17 and July 18 by about 
1 Jy or about 25% of its absolute flux density is remarkable. 
On July 1 8 two observations were made revealing a flux density 
increase by about 0.2 Jy within one hour. At 4.75 GHz and 2.685 
GHz the observed flux densities are similar on July 17 and 18. 
A slightly decreasing flux density at 4.75 GHz within the errors 
is observed. In view of the observational errors given in Table 1 
this decrease is highly significant at more than 10<r at 4.75 GHz. 
Our cross checks on the nearby calibrator 3C138 prove that 
this feature cannot be due to systematic effects. The repeated 
observations on July 18 and July 23 demonstrate the reliability 
of our observations. At 10.55 GHz the source recovered within 
days, at 4.75 GHz within a week, and at 2.695 GHz in less than 
four weeks. 

The same behaviour is seen in the data taken by the Green 
Bank Interferometer. At 8.3 GHz the flux density was extremely 
low at the time we observed the depression at 10.55 GHz with 
the Effelsberg telescope (see also Fig.5). However, since only 
one observation is strongly influenced the feature does not ap- 
pear in the two week averages in Fig. 2b. At 2.25 GHz the flux 
is extremely low for a time of roughly two weeks, and the sig- 



nificance of the depression seen in Fig. 2b is much higher than 
lOcr. This is an independent confirmation that the Effelsberg 
data show a real feature which is only undersampled. The dif- 
ferences in the flux variations at the different frequencies are 
important constraints for the modelling described below. 

The light curves of 0528+134 show that a similar event 
but with a smaller amplitude may have occurred in June 1992. 
There is also some evidence for a counterpart in the Green Bank 
data. 



3. VLBI Observations 

3.1. VLBI observations at 22 GHz, Data analysis and 
Results 

0528+ 1 34 was observed with a global VLBI array of eleven sta- 
tions (see Table 2) at 22 GHz from November 7, 21h30 UT, to 
November 8, 13h30 UT, in 1992. The data were recorded using 
the Mk III VLBI recording system (Rogers etal. , 1983)inmode 
B with an observing band of 28 MHz centred at 22236.99 MHz. 
The correlation was done at the Max-Planck-Institut fur Ra- 
dioastronomie at Bonn. Fringe-fitting and imaging was done 
in the usual way using standard VLBI-software (Alef 1989, 
Pearson, 1991, Shepherd et al. , 1994). To correct for atmo- 
spheric attenuation, the amplitude calibration of the visibilities 
was improved applying opacity corrections to each station (for 
details of the method see Krichbaum et al. , 1993). In Table 

2 the participating stations and their antenna performance are 
summarised. 

In Figure 3 we show the CLEAN-map of 0528+134 obtained at 
22 GHz. The source structure can be described as a one-sided 
bent core -jet structure of ~ 5 mas length. The parameters of 
the individual structure components were determined from fits 
of elliptical Gaussian components to the data. The best results 
were obtained with the 6 component model given in Table 3a. 
In this model distinct jet components were found at separations 
of r = 0.2, 0.6, 0.9, 3.1, and 4.9 mas relative to the brightest 
component. Due to their faintness the two outer components 
at r = 3.1 mas and r = 4.9 mas are less well defined than the 

3 central components and may reflect the possibility that the 
jet emission at separations of r > 3 mas is resolved by the 
interferometer beam and/or just gets too faint (the map contains 
about 90% of the total flux density). During the model fitting 
some uncertainty remained in the determination of the exact 
location of the emission between r = 0.6 mas and r = 1 .0 mas. 
A slightly worse model fit with only 5 components, yields an 
extended component CI /C2 at r = 0.8 mas, which concatenates 
CI and C2 into a single component. In Table 3a the parameters 
of this concatenated component C1/C2 are summarised. The 
non-uniqueness of the models may reflect a relatively sparse 
uv-coverage on intermediate-length baselines. Since the data 
obtained at 8.4 GHz (see below) also showed distinct emission 
near r = 0.5 and r = 1.0 mas, we used the more complex 6 
component model for further discussions. 
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Fig. 1. Lightcurves of 0528+134 from the Effelsberg 100-m and the IRAM 30-m telescopes. The numbers are the corresponding observing 
frequencies in GHz. The errors of the individual data are given in Table 1 . The source activity in the EGRET range is indicated by open 
rectangles for normal level and filled rectangles for a high 7— ray flux. 
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Fig. 2a. (On top) Details of the lightcurves of 0528+134 at 2.695 GHz, 4.75 GHz and 10.55 GHz. Again the 7-ray state is indicated. 

Fig. 2b. (Bottom) The lightcurves of 0528+ 134 at 2.25 GHz and 8.3 GHz as averages over periods of two weeks. For the two weeks in July 1993 
no data at 2.25 GHz are available due to hardware malfunction. The dip at 8.3 GHz in July 93 is averaged out since only one of six observations 
is influenced. 
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Table 1. Flux densities of 0528+134 in 1993 as observed with the Effelsberg 100-m telescope and the IRAM 30-m telescope. 
Frequencies are in GHz and flux densities in Jy. 



Date 
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0.23 
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0.10 


2.62 


0.14 


4.01 


0.20 


7.0 


2.0. 






30.04. 














7.2 


0.5 






01.05. 


1.81 


0.07 


2.57 


0.14 


4.42 


0.30 










10.05. 














6.6 


0.7 






25.05. 














6.1 


0.5 






18.06. 


















8.8 


1.0 


19.06. 


















8.1 


0.9 


20.06. 






2.63 


0.08 


4.29 


0.21 






7.9 


0.9 


22.06. 


















8.7 


1.0 


23.06. 


















7.4 


0.8 


05.07. 
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0.4 
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7.0 


0.4 


17.07. 


1.39 


0.04 


1.68 


0.08 


2.34 


0.16 
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1.31 


0.04 


1.46 


0.06 
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0.15 
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1.29 
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1.45 


0.06 
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0.17 










19.07. 
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0.04 
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4.59 
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0.5 
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0.2 
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2.22 


0.11 


2.87 


0.12 
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0.35 
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7.1 
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12.09. 
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0.3 
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6.9 
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06.11. 






2.50 
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7.8 


1.0 
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0.11 
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0.36 
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0.11 
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0.26 










25.11. 


1.91 


0.07 


2.89 


0.17 


4.48 


0.22 
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4.7 


0.3 
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1.80 


0.10 
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4.72 


0.43 
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0.31 


5.2 
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Table 2. Station characteristics at 22 GHz. Column 1-6 give the 
station name and location, the diameter of the antenna, the typ- 
ical system-temperature during the observations, the peak sen- 
sitivity of the antenna gain, and the mean atmospheric zenith- 
opacity during the experiment, respectively. 



Station 


T oration 


Diam. 


T 

A sys 


(X 
& 








I"ml 


TK1 


TK/Tvl 




LJIClSUClg 




1 on 


1 10 


791 

\J. / 71 


07 


Medicina 


Italy 


32 


150 


0.110 


0.09 


Onsala 


Sweden 


20 


320 


0.053 


0.04 


Green Bank 


West Virginia 


43 


70 


0.141 


0.06 


Hancock 


New Hampshire 


25 


90 


0.087 


0.05 


North Liberty 


Iowa 


25 


100 


0.108 


0.09 


Fort Davis 


Texas 


25 


60 


0.095 


0.06 


Los Alamos 


New Mexico 


25 


90 


0.096 


0.06 


Pie Town 


New Mexico 


25 


100 


0.134 


0.07 


Kitt Peak 


Arizona 


25 


95 


0.113 


0.06 


Owens Valley 


California 


25 


150 


0.108 


0.07 



3.2. VLBI Observations at 8 GHz, Data analysis and Results 

Within the IRIS-S (International Radio Interferometric Survey- 
ing) and EUROPE geodetic VLBI observing campaigns per- 
formed monthly (IRIS-S) and quarterly (EUROPE), the source 
0528+134 is observed regularly in the S- and X-band (2.3 GHz 
& 8.4 GHz). Geodetic VLBI observations therefore provide an 
important astronomical data base to study the mas-structure 
of compact extragalactic radio sources and its kinematics on 
timescales of a few months (Britzen et al. , 1993). From these 
observing campaigns we selected and analysed two data sets 
in the X-band (epochs: February 25, 1991 (1991.23) & March 
25, 1992 (1992.15)), which were relatively close in time to our 
observations at 22 GHz and previous VLBI observations per- 
formed in 1990/1991 by Z94. In these experiments 0528+134 
was observed in snapshot mode for 12hrs using the VLBI an- 
tennas at Mojave (USA), Richmond (USA), Westford (USA), 
Wettzell (Germany) and Harthebesthoek (South Africa). Mainly 
due to the participation of the antenna at South Africa, the 
observing beam was nearly circular with a size of typically 
~ 0.5 mas. The data were correlated and imaged with the same 
standard software as described above. The amplitude calibration 
was done in the standard way, using regular system-temperature 
measurements and antenna gain curves provided by the stations. 
The 'a priori' amplitude calibration based on data provided by 
each observatory was further improved by imaging of structure 
calibrators (sources of simple and known mas-structure, e.g. 
1803+784, 4C39.25, see Britzen et al., 1994) in parallel to 
0528+134, using methods described in Schalinski et al. , 1986. 
The maps and model fits obtained from the two data sets 
for 0528+134 clearly show a one-sided core -jet structure of 
more than 5 mas length, bending from P. A. ~ 60° near 
the core to P. A. ~ 5° — 20° at r > 5 mas. The mas- 
structure at 8 GHz can be best described by 5 distinct features 
CO - C4, which are located (relative to the brightest compo- 



nent CO) at separations of r = 0.34 ± 0.18mas (1991.23)and 
r = 0.57 ± 0.18mas (1992.15), respectively, r ~ 1.1 mas, 
r ~ 3.2 mas, and r ~ 4.9 mas. In Table 3b & 3c we summa- 
rize the model fit parameters of the VLBI components at both 
epochs. 

3.3. Discussion of the source structure 

Inspection of maps and model fits at 8 GHz and 22 GHz reveals 
a core-dominated core-jet structure, with embedded distinct 
features and an overall pronounced and continuous bending of 
about 60° — 70° : the elliptically elongated component CO and 
the innermost jet component (at r = 0.2 mas) are oriented along 
P. A. = 80° , the components near r = 0.6 mas are at P. A. ~ 
60° , the components near r = 1.0 mas are at P. A. ~ 30° , and 
the outer components are located at P. A. = 10° — 25° . 

3.3.1. Identification of the VLBI core: 

The interpretation of the morphology and kinematics of the mas- 
structure of AGN critically depends on the existence of distinct 
jet components, registered at the different observing times and 
frequencies. Little was known up to now on the mas-structure of 
0528+134. A first map published by Chariot (1990) revealed a 
first but only crude estimate of the basic source structure. More 
detailed maps obtained in 1990/1991 at 22 GHz and 8 GHz by 
Z94 yielded two alternative component identifications schemes, 
which basically differed in the interpretation of the brightest 
component as being or not being the VLBI-core of 0528+134. 
With the new data, it is now possible to discriminate between 
these two schemes. 

Based on our results summarised in Table 3 and the maps and 
model fits obtained by Z94 (1990.85: 22 GHz map, 1991.16: 
8 GHz map) we suggest the component identification indicated 
by label 'C in column 1 of Table 3, in which CO is the VLBI 
core, which we assume to be stationary. This identification of 
CO as VLBI core is supported by its pronounced flux density 
variability, its inverted spectrum, the one-sidedness of the mas- 
structure, and the fact that two components (N & CI) seem to 
move away from CO. To include the VLBI components given by 
Z94 in the discussion, we redetermined their positions relative 
to the concatenated components K2K3 and X2X3 (see Table 1 
in Z94), which we identified as the core. Preliminary results ob- 
tained from VLBI observations at 43 GHz and 86 GHz support 
this identification scheme (Standke 1994). 

3.3.2. Motion of CI: 

With the assumption of CO being the stationary VLBI core 
a simple cross-identification of the structure components seen 
in the maps/model fits at 8 & 22 GHz is suggested. In this 
identification scheme the components C2 (at r ~ l.Omas), C3 
(at r ~ 3.2 mas) and C4 (at r ~ 4.9 mas) remained stationary 
relative to CO within the observing interval 1991.16 - 1992.85. 
For component C 1 , however, we find evidence for an increase 
of its core separation with an apparent angular separation rate 
of n = 0.11 ± 0.043 mas/yr (see Figure 4), corresponding to 
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Table 3a. Results of the Gaussian component model fits at 
22.2 GHz, 1992.85. Here and in the tables 3b & 3c column 
1 gives the component identification, column 2 the flux density, 
column 3 and 4 the component position in polar co-ordinates 
(distance, position angle), column 5 the geometric mean size 
of the major axis and minor axis of the Gaussian component. 
Numbers in parentheses are formal errors assuming Gaussian 
error propagation. The position angles are counted from north 
through east. 



Id. 


S 


r 


P. A. 


FWHM 




Vy] 


[mas] 


[°] 


[mas] 


CO 


3.00 (0.05) 








0.08 (0.02) 


N 


0.32 (0.03) 


0.23 (0.03) 


83 (10) 


0.21 (0.10) 


CI 


0.16(0.05) 


0.63 (0.20) 


64 (20) 


0.72 (0.20) 


C2 


0.10(0.05) 


0.92 (0.20) 


30(10) 


0.61 (0.20) 


C3 


0.017 (0.008) 


3.14(0.35) 


27(5) 


0.20 (0.13) 


C4 


0.016(0.011) 


4.9 (0.9) 


23 (10) 


0.76 (0.90) 


Cl/2 


0.22 (0.05) 


0.81 (0.09) 


44(5) 


1.01 (0.19) 



Table 3b. Results of Gaussian component model fits at 8.4 GHz, 
1991.23. 



Id. 


S 


r 


P. A. 


FWHM 




Vy] 


[mas] 


[°] 


[mas] 


CO 


0.79 (0.20) 








0.44(0.10) 


CI 


0.24 (0.09) 


0.34 (0.18) 


58(18) 


0.20 (0.10) 


C2 


0.33(0.15) 


1.13(0.23) 


24(7) 


0.60 (0.20) 


C3 


0.10(0.05) 


3.32 (0.26) 


24(5) 


0.43 (0.30) 


C4 


0.09 (0.06) 


4.8 (0.5) 


6(10) 


0.54 (0.25) 



Table 3c. Results of Gaussian component model fits at 8.4 GHz, 
1992.15. 



Id. 


S 


r 


P. A. 


FWHM 




[Jy] 


[mas] 


[°] 


[mas] 


CO 


1.24(0.07) 








0.37 (0.10) 


CI 


0.21 (0.12) 


0.57 (0.18) 


57(5) 


0.10(0.05) 


C2 


0.10(0.05) 


1.09 (0.10) 


25 (10) 


0.63 (0.20) 


C3 


0.12(0.07) 


3.08 (0.37) 


22(6) 


0.31 (0.15) 


C4 


0.11 (0.07) 


5.0 (0.4) 


12(10) 


0.29 (0.15) 



superluminal motion of j3 app = 4.4± 1 .7 (for z = 2.07 and Hq = 
100 km s _1 Mpc -1 , go = 0.5). Subluminal motion is rejected 
with more than 2 a (4.2% chance probability). The rejection 
is even stronger for a smaller Hubble constant Hq < 100 km 
s -1 Mpc -1 . A linear fit for the angular separation rate has a 
significance level of only 57% in the x 2 -statistics. This is mainly 
caused by the small error bars for the earlier data by Zhang et 
al. (Z94), which partly are the result of the small FWHM of 
the component CI at this time. The error bars provided for 
our VLBI observations were calculated including not only the 
the formal Gaussian errors, but also uncertainties derived from 



0528+134 22.237 GHz 1992.85 




4 2 

Relative R.A. (milliarcsec) 



Fig. 3. Clean-Map of 0528+135 at 22 GHz restored with an elliptical 
beam of size 1.62 x 0.25 mas, pa = — 10° . Contour levels are -0.1, 0.1, 
0.3,0.5, 1,2,5, 10, 15,20,30,50, 70, and 90% of the peak flux density 
of 2.8 1 Jy/beam. The brightest component CO corresponds to the core. 
In its eastern wing the three components N, CI, and C2 are hidden, 
whereas the weaker outer components C3 and C4 can be clearly seen 
in the north-west of the core. The structure east of the core is at the 
noise level. A simple sketch of the inner structure will be shown in 
Fig.6 . 

the internal scatter between slightly different calibrated data 
sets and the resulting slightly different models of the source 
structure, and thus may be rather too conservative. 

Another source of uncertainty arises when using data at 
different frequencies, namely a possible frequency dependence 
of the core separation. For 0528+ 1 34 this effect appears to play 
a minor role as the 8.4 GHz data and the 22 GHz and 43 GHz 
data in Fig.4 follow - within the statistics - the same regression 
line. However, it may be wortwhile to check the stability of 
our results against this systematic uncertainty by enlarging the 
error bars to at least Sr =0.1 mas. Doing so for the two earlier 
data points does not change much of the result, but improves 
the Xmin to 95% significance level. However, then subluminal 
motion is rejected with only 1.8 it (7.3% chance probability). 
A preliminary analysis of 8.4 GHz VLBI observations at epoch 
1993.53 yields an indication that the component CI is a blend 
of two individual components (Britzen et al., in prep.). These 
new observations have been done in the same manner as the 
earlier 8.4 GHz observations except that the stations Santiago 
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in Chile and Fortaleza in Brazil replaced the stations at Mojave 
and Richmond. 

It is noteworthy that the position angle of CI relative to the 
core CO changed smoothly from P.A. ~ 89° in 1990.85 (Z94) 
to P.A. ~ 50° - 60° in 1992 (see Table 3). Since the jet is 
oriented along P.A. ~ 10° — 20° at larger core separations, it 
is tempting to assume that CI follows a bent path. With regard 
to the observed bent jet structure (See Fig. 3) such a behaviour 
would not be completely unexpected and is also seen with mm- 
VLBI in an increasing number of other blazars (e.g. Krichbaum 
et al. , 1994, and references therein). 



1 .o r 



0.8 - 



□ 43.2 GHz 
022.2 GHz 



^8.4 GHz 



£ 0.6 - 



0.4 - 




0.2 - 



0.0 L 



90 



91 



time 



92 

-1900 [yrs] 



94 



Fig. 4. Core separation of component CI relative to CO. Data points 
at epochs 1990.85 (22 GHz) and 1991.16 (8 GHz) are from Zhang 
et al. , 1994. The data point at 43 GHz (epoch 1992.4) is taken from 
Krichbaum et al. (1 995). We have further added a data point from a pre- 
liminary analysis of 8.4 GHz observations at epoch 1993.53 (Britzen 
et al, in preparation). The line represents a x 2 -fit to the data with 
parametrization r(t) = (0.11 rt 0.043) • (t - 1990) + (0.32 ± 0.05). 



3.3.3. Ejection of a new jet component ? 

The 22 GHz data revealed the existence of a new component at 
r = 0.23 ± 0.03 mas, east of the core CO (P.A. = 83 ± 10° ), 
not seen previously at 8 GHz or 22 GHz, but close to the 1990- 
position of its precursor CI. We have labelled this component 
'N' (see Table 3a). This component is also seen in a 43 GHz 
map from 1992.40 (Krichbaum et al. , 1995) and at 86 GHz 
from 1993.3 (Standke, 1994). It is therefore likely that a new 
jet component was ejected prior to 1992.4 and is now moving 
outwards. Disregarding possible small frequency shifts due to 
opacity effects between 22, 43, and 86 GHz, we tentatively 
determined the angular separation rate of 'N' to be about < 
0.5mas/yr (f3 app < 20). Future observations should clarify, 
whether component 'N' was ejected during the high-frequency 
radio outburst in 1992.1. 



3.3.4. Spectral indices of the jet components: 

With the available data it is possible to determine the spec- 
tral shape of the jet components. In Table 4 we summarise 
the results for the spectral indices a%/22GHz an d "22/43 GHz- 
The errors for the spectral indices were calculated formally, 
using the errors of the corresponding flux densities given in 
Table 3 and assuming Gaussian error propagation. Between 
1990/1991 and 1992 the spectrum of the core component CO 
inverted considerably, reflecting that the flux density of the core 
increased - during an ongoing outburst - first at the higher fre- 
quency. With spectral slopes of a%/22GHz = —0.90 it 0.10 and 
a 22/43GHz = 0.30±0. 10a spectral turnovernear v m = 22 GHz 
could be expected for CO. However, the variability during 1992 
in our 32 GHz data and the 22 GHz and 37 GHz data in Z94 
imply that this determination of a turnover frequency should 
not be taken too seriously. 

The apparent spectral steepening of the component spectra 
with increasing core-separation (see Table 4) has to be regarded 
with some care. Whereas the spectral slopes of the inner jet 
components N, CI, and C2 are smaller than 0.7, the spectral 
indices of the components C3 and C4 lie in the range of 1 — 2. 
At 22 GHz, however, the outer components C3 & C4 are most 
probably resolved by the interferometer beam, leading to under- 
estimates of their flux density and therefore to upper limits of 
the corresponding spectral index a 8 / 2 2- Within still relatively 
large uncertainties the spectral slopes for the jet components 
N-C4 therefore are consistent with typical optically thin syn- 
chrotron spectra (as observed in comparable objects), and gives 
us confidence that our interpretation and cross-identification of 
the features seen at 8 GHz, 22 GHz, and 43 GHz is not implau- 
sible. 



4. Discussion 

4.1. Luminosity limits and superluminal motion 

An interesting aspect of the analysis of AGN 7-ray data is 
to look for luminosity limits above or below which isotropic 
emission from the source is not possible. A violation of these 
limits can then be taken as argument in favor of an origin of 
the emission in a relativistic jet pointing nearly towards us. In 
such a scenario superluminal motion of VLBI components is 
expected. 

Under the synchrotron-self-Compton model Zhang et al. 
(Z94) have tried to determine the physical conditions of the 
emission region on the basis of ROSAT data taken in March 
1991. If the angular size of the emission region is known (from 
VLBI) and one is able to get the turnover frequency between 
optically thick and optically thin synchrotron emission from 
fits of the radio spectrum, then the self-Compton X-ray flux can 
be calculated. Z94 found that the predicted X-ray flux is more 
than three orders of magnitude higher than the flux observed 
by ROSAT, and conclude that Doppler boosting is required in 
0528+134. These authors used a turnover frequency of 7 GHz 
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Table 4. Spectral indices of the components. Column 1: component identification, column 2: observing epochs of fre- 
quency 1 /frequency 2 used for aj requency i/ j requency 2, column 3: core separation of the component at epoch 1992.85, column 4 
& 5: spectral indices a%/22GHz an d a 22/43GHz, using S oc v~ a . 



Id. 


interval 


r(1992.85) 
[mas] 


a i/21GHz 


"22/43 GHz 


comment 


CO 


1990.85/1991.16 





0.30±0.10 




from Z94 


CO 


1992.15/1992.85 





-0.90 ±0.10 






CO 


1992.85/1992.40 







0.30±0.10 


using 43 GHz data 


N 


1992.40/1992.85 


0.2 




0.70 ±0.70 


using 43 GHz data 


CI 


1992.15/1992.85 


0.6 


0.30 ±0.70 






C2 


1992.15/1992.85 


0.9 


0.10±0.10 






C3 


1992.15/1992.85 


3.1 


1.6 ±0.5 




upper limit, see text 


C4 


1992.15/1992.85 


4.9 


1.9 ±0.9 




upper limit, see text 



and a VLBI size for the core of 0.17 mas, which they obtained 
from 22 GHz VLBI observations taken half a year earlier. Both 
parameters are very crucial since they enter the prediction by 
the sixth and eighth power respectively. We have reinvestigated 
this case and find that at the time of the ROS AT observations the 
VLBI size at 8.4 GHz was 0.4 mas, both in our observations and 
in observations by Zhang et al. (Z94), which they did not use for 
their analysis. This size should be taken as an upper limit, since 
the core is not resolved; however, it is the only information we 
have from that time. 

Furthermore, at the time of ROSAT observations the radio 
flux from 0528+134 was around 2.1 Jy at 2.25 GHz, 1.7 Jy 
at 8.3 GHz (our NRL-GBI data) and 2.4 Jy at 90 GHz (Reich 
et al. 1993). Other at least nearly contemporaneous data are 
not available. Therefore, the turnover frequency in the radio 
spectrum can not be easily determined. It is, however, surely 
not 7 GHz. One could argue for a optically thick subcomponent 
in the core. But also a turnover frequency of 80 GHz with a 
size of 0.1 mas give a prediction which matches the ROSAT 
flux. Thus Doppler boosting is not required in view of the low 
luminosity in the ROSAT range. The assertion by Z94 is no 
longer tenable. 

Recently, McNaron-Brown et al. (1995) argued for beaming 
of the MeV emission from 0528+ 1 34 on the basis of OSSE data. 
The luminosity and the variability timescale at these energies 
implied a violation of the compactness limit, i.e. the optical 
depth against photon-photon pair production was much larger 
than unity. 

Due to the strong reduction of the Compton cross section in 
the Klein-Nishina limit the Eddington limit in terms of the lumi- 
nosity and the variability timescale (Elliot and Shapiro 1974) do 
not force us to the conclusion that the GeV-emission of AGNs is 
strongly Doppler-boosted. Also the pair production cross sec- 
tion is reduced at photon energies in the EGRET range and 
hence the corresponding optical depth is less than unity for the 



interaction between high energy 7-rays (see also Dermer and 
Gehrels 1995). However, the compactness limit was violated 
by 0528+134 in the medium energy 7-ray range (1-30 MeV), 
thereby indicating a requirement for Doppler-boosting. A de- 
tailed discussion of luminosity limits in the relativistic range, 
i.e. for 7-ray energies exceeding the rest mass of electrons, can 
be found in Appendix B. 

Our finding of superluminal motion nicely complements to 
the requirements for relativistic motion due to the luminosity 
limits. Therefore it provides strong evidence that blazar radia- 
tion is produced by a relativistically outflowing jet which beams 
emission towards favourably oriented observers. 

4.2. The lightcurves 

In March 1993 0528+134 exhibited an extremely strong 7— ray 
outburst (I(E > 100 MeV) ~ 3 • 10" 6 cm" 2 sec _1 ) lasting for 
a couple of days (Michelson et al. 1994). In May 1993 the 7-ray 
flux was down to a normal level similar to the flux observed in 
fall 1992. We noticed no simultaneous radio flare at that time 
except for a secondary maximum at 32 GHz. On the basis of 
monitoring data at 22 Ghz and 37 GHz Z94 concluded that the 7- 
ray flare occured just prior to a sharp peak at radio frequencies. 
With our data we can say that this maximum in the radio flux 
density was a secondary one and that it occured between the 
two 7-ray observations, the first showing an outburst and the 
second showing low-level emission. In fact the 32 GHz flux 
density is the same within the errors at the exact time of EGRET 
observations (see Table 1 ). The 22 GHz and 37 GHz data by Z94 
end roughly at the time of the main radio flare when on average 
the flux density is higher than for the secondary maximum in 
April except for one data point at each frequency which may be 
either flickering or ill-determined. 

It may well be that a simultaneous flare occured at optical or 
NIR wavelengths, but we prefer to draw a relation between the 
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extremely bright 7-ray flare and the main radio flare which be- 
came visible a few months later during the summer of 1993 and 
not between the 7-ray flare and the small secondary maximum 
in the radio lightcurve. 

We have speculated that the new VLBI component N could 
have been expelled from the core during the outburst at high 
radio frequencies in the beginning of 1992. Since in summer 
1993 the intense radio flare at all frequencies above a few GHz 
appears to be much stronger than the outburst in the beginning 
of 1992 it may well be that again a new VLBI component was 
born, which could be checked by new VLBI observations. 

We found a rather unusual strong depression in the radio 
light curves in July 1993 when the source flared at 32 GHz 
and 86 GHz. Based on the data presented above we discuss two 
cases as the most likely explanation for the this strong variability 
of 0528+134. Firstly we consider a possible origin within the 
source itself. A different approach suggested by the permanent 
fluctuations at 2.25 GHz and 2.695 GHz and the dark cloud 
Barnard 30 along the line of sight is the effect of an extreme 
scattering event by the interstellar medium. 

4.3. Intrinsic variations 

Obviously, any modelling of the July 1993 event suffers from 
missing data coverage in the first half of the month. We can 
therefore not exclude the possibility that the event is caused 
by intrinsic properties of the source itself. The strong 7-ray 
emission and variation on very short time scale as well as the 
available VLBI observations indicate that a relativistic jet is 
pointed towards us with a very small deviation from the line 
of sight. In such a geometry even very small bending effects 
may cause large amplitude fluctuations as has been studied for 
the case of an anisotropic jet with non-spherical emission re- 
gions by Reynolds (1982) and Reynolds and Ellison (1992). 
Such effects may even lead to intraday variability in the optical 
synchrotron radiation and in the radio regime. However, since 
this process could influence only one VLBI component at a 
time, it appears difficult to reproduce the observed behaviour 
by changing the viewing angle towards the flaring core com- 
ponent CO. Judged by their intensity in the 22 GHz VLBI map, 
all components except the core do not provide enough flux at 
10.55 GHz to reduce the total radio emission by 2 Jy if they are 
switched off. Different components should not know what the 
others do and thus intrinsic variations have to be related to the 
core component. Due to its optical thickness it is questionable 
whether the core also provides enough 2.25 GHz and 2.695 GHz 
flux to explain the observed behaviour. However, at 32GHz and 
86GHz no variation is observed within an uncertainty of around 
10% which can hardly be accounted for in emission models of 
an anisotropic jet. Additionally, component CO would have to 
return to its original path after the event, which is an unlikely 
behaviour for a bent jet. 

One could also think about an optically thick synchrotron 
blob in the jet passing in front of the core and thereby obscuring 
part of its emission. The temporal behaviour of the event would 
then have something to do with frequency dependent substruc- 



ture in the core. However, a transverse motion by 0. 1 mas in 
about a week implies velocities of around 200 times the speed 
of light which is impossible. Please note that here the effect of 
apparent superluminality does not occur. 

A third possibility was discussed by Marscher (1979) on the 
basis of an unpublished suggestion by McAdam. The basic idea 
is that a layer of HI in front of the source gets ionized by UV 
photons from an optical outburst and then provides an optically 
thick screen due to free-free absorption. However, covering the 
core and the inner components with such a screen requires fine- 
tuning to the 10~ 4 level, which is the duration of the event 
divided by the light travel time for the system. Furthermore, the 
required emission measure EM~ 10 8 pc cm -6 and an extremely 
dense layer is required to get the recombination time scale to 
the order of a day. 

To summarise, we regard the possibility of an intrinsic cause 
for the observed variability as unlikely. 

4.4. Extreme scattering in Barnard 30 

In the last decade several studies have focused attention on re- 
fractive and diffractive effects of electron density fluctuations 
in the interstellar medium upon radio observations of pulsars 
and compact extragalactic radio sources. Both low-frequency 
intensity variations and time-of-arrival fluctuations in case of 
pulsars have been successfully modelled as scintillation due to 
interstellar turbulence (Sieber 1982; Rickett 1986; Cordes, Pid- 
werbetsky and Lovelace 1986; Romani, Narayan and Blandford 
1986; Goodman et al. 1987; Rickett et al. 1994). While these 
effects are stochastic in nature, the so-called extreme scatter- 
ing events have provided evidence for the existence of discrete 
large-scale inhomogenuities in the ionised interstellar medium 
(Fiedler et al. 1987a; Cognardet al. 1993). 

We have already noted that the variability of about 25% at 
2.25 GHz and 2.695 GHz on the timescale of a few months gives 
some evidence for stochastic interstellar scattering (RISS). It is 
therefore not unlikely that from time to time extreme scattering 
events also occur. 

The mathematical details of our model can be found in Ap- 
pendix A. The physics involved here is basically the dispersion 
of electromagnetic radiation in a plasma. The phase velocity 
of radio waves in a plasma is no longer equal to c, but slightly 
larger. The scatterer thus imposes a phase shift on the traversing 
waves similar to a diverging lens. The resultant intensity profile 
can then be calculated following Huygens' principle, that each 
element of the distorted wavefront is the center of a secondary 
disturbance which gives rise to spherical wavelets. The position 
of the wavefront at any later time is determined as the envelope 
of all such wavelets. 

The model curve in Fig. 5 is calculated for a cloud of ionised 
material with a soft density profile, such that the line-of-sight 
integral of an electromagnetic wave passing at distance p from 
the cloud centre follows 

/ al n e = — 

J_oo r + p 
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Fig. 5. Observed and modelled light curves at 2.25 GHz (crosses, 
scaled down by 0.3), 2.695 GHz (squares, scaled down by 0.5), 4.75 
GHz (plus signs, scaled down by 0.8), 8.3 GHz (asterisks), and 10.55 
GHz (diamonds) in July 1993. Source evolution during the event is 
taken as a smooth increase by about 20% at 8.3 GHz and 10.55 GHz 
peaking on the 21st of July, in order to agree with the low flux level at 
8.3 GHz in the beginning of July 1993. The model light curves for an 
extreme scattering event harmonise well with the data. 
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Fig. 6. This plot shows the path of the occulter relative to the inner four 
VLBI components of 0528+134 as used in Fig.5. The components N 
and C 1 have been placed corresponding to their apparent velocity and 
the time lag of 0.7 years between the 22GHz VLBI observations and 
the scattering event. The circle's sizes represent the FWHM-sizes of 
the components, which in our modelling were taken to be independent 
of frequency. The size of the circle drawn for the occulter represents 
the size scale ro and not the overall size of the scattering system. 



where 

n e =400crrr 3 , r = 4-10 11 cm. (1) 

These numbers are similar to those which Cognard et al. (1993) 
have deduced for an extreme scattering event in the direction of 
the pulsar 1937+21, although they used a much harder density 
profile. Due to the time consuming nature of the numerical 
analysis of the integral (A13) these numbers are not the result 
of a fitting process. The fits and the parameters we give therefore 
illustrate reasonable, not optimal, fits to the data. 

We assume this cloud to be a local condensation in the 
environment of the dark cloud Barnard 30, thereby taking its 
distance to be 400pc. It has been noted before that extreme 
scattering events occur preferentially near galactic foreground 
structure (Fiedler et al. 1994). 

The path of the scattering cloud relative to the inner VLBI 
components of 0528+134 is shown in Fig.6. The outer com- 
ponents C3 and C4 are outside the range of this plot, but still 
included in the model. The components N and CI have been 
shifted according to their apparent motion and the time lag of 
0.7 years between our 22 GHz VLBI observation and the event 
discussed here. The size of all components is taken as indepen- 
dent of frequency and equal to 0.1 mas. A larger size of the 
components would result in a slight smearing of the modelled 
light curves. The apparent transverse velocity of the occulter 
is around 70 km/sec, of which part is due to the motion of the 
earth around the sun. 
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Fig. 7. Here we show the spectra of the VLBI components which we 
used for our modelling. The only significant difference to the observed 
flux densities of the VLBI components is that the core CO is assumed 
to undergo a major outburst in July 1993. 

In Fig. 7 we plot the spectra of the VLBI components as 
taken for our modelling. We assume no significant spectral 
changes compared to the VLBI results at 22 GHz and 8.4 GHz 
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except for the flaring core component CO. The spectra of all 
components are taken as stationary for our modelling. For com- 
ponent C2 we are not sure what the spectrum is. While being 
rather strong at 8.4 GHz in our VLBI data from 1991.23, it does 
not show up in data at the same frequency taken a month earlier 
(Z94). The spectrum we use here fits into the time evolution 
scheme of a cooling electron spectrum in the components. The 
total flux of the source is chosen to match the observed flux 
values before (June) and after (August) the event. To do so we 
have to assume that source evolution increases the 8.4 GHz and 
10.55 GHz flux densities smoothly (with a rise time of 10 days) 
by roughly 20% peaking around the 21st of July following the 
general trend visible in Fig.2a and 2b. 

The observed flux variations in July 1993 can be nicely 
reproduced by our model. This is especially true for the time 
lag between the minimum at 10.55 GHz (before 18th of July) 
and the minima at lower frequencies (after 18th of July). This 
behaviour is not due to our choice of the density profile in the 
scattering cloud, but very much related to the VLBI structure 
of the source. While passing near the line-of-sight the scat- 
tering plasma cloud first influences the emission with inverted 
spectrum from the core component which leads to the early 
minimum at 10.55 GHz. A few days later the emission from the 
inner jet components, which have flatter spectra, is scattered 
yielding later minima at low frequencies. Small fluctuations in 
the density profile of the scatterer show up preferentially at 
2.25 GHz and 2.695 GHz and may cause deviations between 
the model and the real data, whereas at higher frequencies only 
the general density profile and the VLBI structure of the source 
play a role. 

We regard the observed dips in the light curves of 0528+1 34 
as likely to be explained by an extreme scattering event. It may 
well be that similar events in the direction of 0528+134 are 
rather common due to turbulent substructure in Barnard 30. 



5. Conclusion 

In this paper we have reported on multifrequency radio obser- 
vations between 1992 and 1994 of the 7-ray quasar 0528+134. 
In summer 1993 we observed a major outburst at all frequen- 
cies higher than a few GHz, which peaked a few months after a 
strong outburst in high energy 7-rays and showed the canonical 
time evolution of a cooling and/or expanding electron distribu- 
tion. Our data indicate that the outburst in the EGRET range 
originated very close to the central object of the AGN and that 
a remnant of this outburst moved further outward in the jet until 
it became optically thin at radio frequencies after a few months. 

We have also presented a new VLBI map at 22 GHz of 
0528+134 observed in November 1992. The source structure 
is similar to previous VLBI maps except for a new component 
which is not included in a 22 GHz map taken in November 1990 
(Z94) but is seen in 43 GHz data from May 1992 (Krichbaum 
et al. 1995) and at 86 GHz in April 1993 (Standke 1994, PhD 
Thesis). The apparent transverse velocity of this new component 
is Pa VV <20. 



A second component can also be identified in the 22 GHz 
data from 1990 (Z94) and in 8.4 GHz maps. Its apparent trans- 
verse velocity is /3 app =4.4±1.7 (Ho=100, qo=0.5). This result 
gives evidence for the existence of superluminal motion in the 
EGRET source, which was expected since Doppler boosting is 
required to satisfy the compactness limit in terms of medium 
energy 7-ray luminosity and variability timescale. Since only 
few EGRET sources have been looked at with sufficient reso- 
lution in VLBI, the high percentage of superluminal sources in 
the sample of EGRET detected quasars yields strong evidence 
for the assertion that high energy 7-ray emission from AGN 
originates in relativistic jets and is strongly boosted around the 
direction of the jet. 

In July 1993 we observed an unusual strong decrease of 
the flux density by about 50% at 4.75 GHz and 10.55 GHz and 
slightly less at 2.695 GHz wavelength. This feature can also 
be found in the monitoring data at 2.25 GHz and at 8.3 GHz 
from the Green Bank Interferometer and must be taken as real. 
An intrinsic origin of the variability behaviour is found to be 
unlikely on the basis of the source structure, the spectra of its 
components and the frequency dependence of the dips in the 
light curves. 

Alternatively, we model the event by extreme scattering on 
a plasma cloud at 400pc distance near Barnard 30. Extreme 
scattering events are often found in pulsar light curves in which 
the time-of-arrival variations allow a good distinction between 
scattering and other effects. We have complemented our light 
curves by data at 2.25 GHz and at 8.3 GHz from the NRL-GBI 
program which is designed to investigate extreme scattering 
events. It is found that the structure of the source is more im- 
portant for the model light curve than the shape of the scattering 
plasma cloud is, which in our model was described by a soft, 
cylinder-symmetric density profile. 

While for the stochastic interstellar scintillations the impor- 
tance of the VLBI structure of the background source for the 
resulting light curves has already been pointed out by Spangler 
et al. ( 1 993), in this paper for the first time a model of an extreme 
scattering event is based on the observed source structure and 
the spectra of the VLBI components of an AGN. It therefore 
provides the best available evidence for the existence of extreme 
scattering events in the light curves of AGN and may together 
with the statistical information from the NRL-GBI monitoring 
program help to give insights into the sub-AU scale structure of 
the interstellar medium in our Galaxy. 

Appendix A: Radiation transport through ionised matter 

The phase shift of an electromagnetic wave due to the passage 
through a cloud of ionised matter at distance p to the cloud 
centre is 

HP)= 2 -^ dln e = ^- , (Al) 
" J -00 r o + P 

where we have assumed an axisymmetric density profile 

/ al n e = — 

J_oo r + p 
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In reality such a profile can not be regarded as extending to 
infinity, but will be immersed in small scale density fluctuations. 
The density profile (Al) is still a good approximation if it 
dominates the true profile including fluctuations only at radii 
smaller than p d where 



p d > (2Trn e \ 2 zr e rl) 



2\ 1/3 



3 • 10 



12 



A 



2/3 



llcm J 



cm (A2) 



for our parameters. Our description will of course also not work 
for x > pd, i.e. when the line-of-sight to the AGN traverses 
regions far from the scatterer, which are dominated by small 
scale density fluctuations, around a month before or after the 
event. The mass contained by the scatterer is 



M, 



7 ■ 10' 

fion 



^ - In f 1 + ^ 



grams 



043) 



Here / son denotes the degree of ionisation. 

With the wave propagator in the Fresnel regime (Born and 
Wolf 1975) we obtain for the amplitude function (relaxing the 
normalisation) 



A(x) = 



exp 



l\z 



ds 



s j. $ i 

COScp -=. 

s r . rj+^/s 



(AA) 



where x is the distance at the cloud between the observer's 
line-of-sight to the AGN and the symmetry axis through the 
AGN and the cloud, z is the distance between observer and the 
scatterer, 8 = zg,H> 1 = ^p> , r f = r ^ and s c = ^ . 

We use the method of stationary phase for the s-integral 
(e.g. Wong 1989). 



f(s) = s- J— cos 

f'(s) = 1 



Sc 

cose 



$1 



2^/scli 2^/s(r-f + ^/s) 2 



(A5a) 
(A5b) 



/"(*) = 



cos<j> 



-3/2 



$1 



$1 



4s 3 / 2 (r; + ^/s) 2 2s(r s + ^/sf 



045c) 

The condition f'(s) = gives a polynomial of third order in 

y = \fs- 



y s + 2r f 



COS0 

2\fs~ c ~ 



$j rjcos<i 



2\fs~ c ~ 



7 I , rt cose 
y z +[rj- 



■ 



2r, 



With the variable x = y + ,. 
equation 

x 3 + px + q = 



(^46) 



■tr^L we obtain the reduced 



(Al) 



where 



X COS<7 



- — r) [l + 3 £ + 3e 2 + e 3 



= 27$o 

27 1 L /J ' ArQr 2 

For the discriminant D = (p/3) 3 + (q/2) 2 we derive 
D = 3~ 6 r 6 s [q 2 -2r 1 {l + ef] . 



(AS) 
(A9) 

0410) 



Depending on the sign of D we have one or three real zeros y 
of the above equation. Interesting for us are only the positive 
ones. The ^-integral is then carried out only for those $ for 
which a positive zero exists. 

With the reduced zero j/q = fj we get for the coefficients 
of the expansion 



/(so = 2/ ) = r ] \Vo 
and 



2e Vo 



/"(so) = r 



S 



27 1 + y* 
7/0+M) 



(All) 



2yf 27 j/q 3 (1 + 2/q) 3 



= r f (7(2/0) • 
For the amplitude function we then get 



0412) 



A(x) ~ V2r / rf$ 
Jo 



V— expfi[/($)±-1 



(A 13) 



which can be solved numerically. Here the sign is determined 
by sgn#($) in the sense that the argument of the square root 
must be positive. 

Appendix B: Luminosity limits in the extreme relativistic 
range 

There are basically two luminosity limits for AGN, the Elliot- 
Shapiro relation and the compactness limit. The former is the 
Eddington limit written in terms of the variability timescale r 
as light travel time for the Schwarzschild radius of the source. 
In the relativistic limit the Thomson cross section is no longer 
valid and one has to use the Klein-Nishina cross section. Then 



L < 



a T ln(2e) 



(51) 



This limit is not violated during the outburst of 0528+134 in 
March 1993, for which an integrated photon flux of I(E > 
100 MeV) ~ 3 • 10" 6 cm" 2 sec" 1 is reported (Michelson et al. 
1994). The spectral index was s = 2.25 ±0.11 and there is 
weak evidence for short-time variability. This limit was also 
not violated in May 1991 when the source was detected for the 
first time. At that time the integrated photon flux was a factor 
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2.5 less, the spectral index was about 2.6, and clear evidence 
was found for variability on timescales of less than a day in the 
source frame (Hunter et al. 1993). 

Now we discuss the compactness limit in the relativistic 
range. The optical depth of photons with energy c\ (here all 
energies are in units m e c 2 ) due to pair production by photon- 
photon interactions in a spherical emission region of radius R 
and differential photon luminosity L(e) is 



/(ei) ~Rf dtt(l - (j) I de' n'(e', ft) a(e') 
where 



n'(e',fi): 



4e' 



L 



16tt 2 R 2 c ei (1 - n) V £ i (!"/") 



2e> 



a 



(52) 



(53) 



is the differential photon density in the CMS of the photon- 
photon system and 



.( £ ') = ^(l-/3 2 ) 



with 



(3-/3 4 )ln|^-2/3(2-/3 2 ) 



(54) 



is the total cross section for pair production in the CMS (Jauch 
and Rohrlich 1976). 
In case of a spectrum 



L — Lq E m in £ , 

we get 

3io5 ms ' n ei 



e > E„ 



8tttc 2 



d,i (i - 



de 



(55) 



(56) 



with 



„ = max (1.0, \/0.5E min e x {\ - (j) 



and r as the minimal variability timescale. In a go = 0.5 cos- 
mology we can write the luminosity in terms of the observed 
differential photon spectrum 1(E) = I E' min E~ 2 



LqE„ 



16^c 2 (C" a - 1) 2 Jq5„ 
Hi 



(57) 



where 

C = 1 + z ; E m i n = (,E min 

Then 

6(C 05 - \fhE min t x 



1(a) 
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Fig. 8. Compactness of an AGN in dependence of the low energy cut- 
off of its gamma ray spectrum. Here E mtn is in units of m e c 2 . The 
compactness scales linearly with the _E_F(_E)-luminosity in units of 
10 49 ergs/sec and inversely with the variability timescale in the frame 
of the AGN in units of days. For 0528+134 the scale factor is around 
unity when Ho = 100 km s _I Mpc _I , go = 0.5 is assumed. 

In Fig. 8 we have plotted the integral (B8) for different low- 
energy cut-offs in the emission spectrum of 0528+134. There 
are 2.5<r detections at 3-30 MeV from COMPTEL for the view- 
ing periods in May 1991 and in March 1993 (Collmar et al. 
1994 and private communication), which indicate a continu- 
ation of the EGRET spectra down to a few MeV and flatter 
spectra at lower energies. In May 1993 0528+134 was not de- 
tected by COMPTEL. One should note that when adding all 
data 0528+134 is detected by COMPTEL with a significance of 
more than 5<t, so it is clearly there, although for single viewing 
periods there are three periods with ~2.5<t detections and a few 
with less than la. It is thus reasonably probable that the three 
~2.5<t detections are real, although the significances for the sin- 
gle observations are very small. The nonlinear dependence of 
the significance of a source on the square root of the observing 
time is due to COMPTEL-specific background problems. 

One should note that on the basis of EGRET data alone one 
would not be able to see the signature of such an attenuation 
since in an optically thick situation the effective optical depth 
depends on the attenuated spectrum of photons. The cut-off 
energy can be roughly determined from Fig. 8. It is that value 
of E m i n for which the compactness at photon energy 5 msn is 
/ = 1. In our case this is 5 msn c ~ 50 which corresponds to 
E c ~ 8 MeV in the observer frame. 

Since the variation of the photon flux in the COMPTEL 
range follows the general trend seen by EGRET, it is natural 
to assume that the 7-rays in both energy regime originate in 
the same region. From Fig. 8 we can deduce that the photons 
with energies less than 10 MeV (observer frame) should have 
been strongly attenuated if the emitting regions were station- 
ary, which is in contradiction to their detection by COMPTEL. 
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Therefore, Doppler boosting of the medium energy 7-ray emis- 
sion of 0528+134 is required. 
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